arXiv: 1509.02393v2 [cond-mat.stat-mech] 24 Sep 2015 


Role of trapping and crowding as sources of negative differential mobility 


Marco Baiesi,^’Attilio L. Stella/’^ and Carlo Vanderzande^A 

^Department of Physics and Astronomy, University of Padova, Via Marzolo 8, 1-35131 Padova, Italy 
‘^INFN, Sezione di Padova, Via Marzolo 8, 1-35131 Padova, Italy 
^ Faculteit Wetenschappen, Hasselt University, 3590 Diepenheek, Belgium 
*Instituut theoretische fysica, Katholieke Universiteit Leuven, 3001 Leuven, Belgium 

(Dated: September 25, 2015) 

Increasing the crowding in an environment does not necessarily trigger negative differential mo¬ 
bility of strongly pushed particles. Moreover, the choice of the model, in particular the kind of 
microscopic jump rates, may be very relevant in determining the mobility. We support these points 
via simple examples and we therefore address recent claims saying that crowding in an environment 
is likely to promote negative differential mobility. Trapping of tagged particles enhanced by increas¬ 
ing the force remains the mechanism determining a drift velocity not monotonous in the driving 
force. 

PACS numbers: 05.40.-a, 05.70.Ln 


I. INTRODUCTION. 

In a system at thermodynamic equilibrium, a particle 
responds to a gentle push by moving on average toward 
the direction of the applied force. Therefore the mean- 
velocity to force ratio, called mobility, is a positive quan¬ 
tity. In general, linear response coefficients are positive 
in equilibrium because the structure of the fluctuation 
response relations involve positive correlations between 
observables and the entropy that perturbations would 
produce [T]. Far from equilibrium there is more vari¬ 
ability. For example, thermal conductivities that are not 
monotonous in the external forcing have been found [2]. 

If a pushed particle starts to go slower by increasing 
the force we have the phenomenon of negative differen¬ 
tial mobility (NDM) [SHU] If the particle responds by 
drifting against the force, we speak of absolute negative 
mobility. The latter behavior has been observed in mod¬ 
els subject to ratchet effects m, in systems driven by 
periodic forces |14[ I15j and for driven Janus particles in 
corrugated channels HU. 

Here we focus on NDM. Negative differential conduc¬ 
tivity was discovered in the 60’s and realized in Gunn 
diodes and other semiconductors at low temper¬ 
ature HIHIS]- Other examples of NDM are listed in 
Ref. [in]. A recent application of NDM is the sorting 
of soft matter colloidal particles [20] . 

For common overdamped diffusion, NDM is expected 
when particles’ surrounding environment contains shal¬ 
low barriers, called “traps” hereafter, in which the par¬ 
ticle driven by an external force might spend more and 
more time upon an increase of the force. For this reason, 
the mobility of the particle may actually decrease rather 
than increase if the force is raised. An example of this ef¬ 
fect was given by particles diffusing within a percolating 
cluster [3] (an instance of hopping process in disordered 
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media EDUl]), in which traps are represented by the 
dangling ends and all portions with concavity against the 
direction of the force. A pedagogical exemplification of 
this “getting more from pushing less” phenomenon was 
then given by Zia and coworkers |4]. In a two-lanes jump 
system that can be solved exactly, NDM emerges when 
one lane contains a hook that can trap the particle. The 
NDM may be described liiniiii] from the point of view 
of a recently developed linear response theory for general 
nonequilibrium systems [33It25j . In this context, NDM 
emerging in simulations of ring polymers impaled by rods 
of a lattice with defects m might explain irregular mi¬ 
gration speeds in experiments comparing linearized and 
circular DNA plasmids HSIIIT]. NDM was also found for 
particles diffusing in a crowded environment mini HO]. 

One expects to observe NDM if there is a form of 
trapping due to a specific coupling between the tracers 
(tagged particle, polymer, etc.) and their environment. 
Jamming of particles, such as in kinetically constrained 
models inis], may also be seen as a form of mutual trap¬ 
ping. However, just a density of obstacles does not say 
anything straightforward about such propensity to meet 
traps if the features of the obstacles and of the tracers are 
not explicitly specified. Yet, using the formalism of meso¬ 
scopic jump systems, it was recently proposed m that 
NDM should be a generic feature for particles strongly 
pulled in crowded environments. 

We believe that shifting the focus from the trapping 
mechanism to the role of crowding is not functional to 
a better understanding of NDM. To support our point, 
in this paper we show that some of the results presented 
in [To], and more in general the appearance of NDM, 
may depend crucially on the choice of microscopic jump 
rates. Moreover, we use some simple models to show that 
crowding may be uncorrelated with NDM or even anti¬ 
correlated. This occurs when too many tagged particles 
are introduced in the system, so that traps are saturated. 
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II. JUMP RATES AND A BASIC EXAMPLE. 

As in previous studies, the results are exposed in the 
context of Markov jump processes. The states of the sys¬ 
tem are discrete and evolve with jumps taking place in 
continuous time. More specifically, we have in mind sub¬ 
sets of the square lattice where the particle may occupy 
free sites and jump to the first neighbors if the transitions 
are allowed. Forbidden transition may occur if the target 
of the jump is already occupied or if there is an idealized 
wall separating the two sites involved in the jump. We 
use a force F parallel to the cc-axis and periodic bound¬ 
ary conditions in this direction. Hence, we monitor the 
average velocity v = {v)p in the x direction. The dif¬ 
ferential mobility ^ = dv/dF becomes negative in the 
regions where v decreases for increasing F. 

For the sake of simplicity we use a unit temperature 
times Boltzmann constant combination, ksT = 1, and 
unit spacing between sites of the square lattice. This 
means that the principle of local detailed balance |2H] is 
met if the jump rate k{C —>■ C) is related to the jump 
rate k(C' —>■ C) of the reversed transition by 

f|^=»pP(C)-C/(C)| (1) 

where U{C) is the energy associated to the configuration 
C. When a particle is subject to a constant force F and 
the jump C —?► C' is in the direction of the force, the 
energy difference becomes the work done by the force, 
U{C) — U{C') = F. When the jump is orthogonal to F, 
simply U(C) = U{C'). There are many (actually infinite) 
choices of the jump rates that satisfy the local detailed 
balance condition. 

Normalizing jump rates of all four exit directions from 
a site with a common denominator is a popular choice 
(see e.g. mziiin]), possibly because it is elegant and re¬ 
sembles a partition function normalization of Boltzmann 
weights. As a consequence, the rate of jumps directed 
transversely to a large force F scale as ~ exp(—F). An¬ 
other possibility would be to keep transversal rates inde¬ 
pendent of F [5]. We call model A the latter and model 
B the former. The details of these rates are given in Ta¬ 
ble |lj where we also specify the scaling for large F as a 
function of e = e~^l'^. We pick up just these two variants 
among infinite ones because they already represent two 
distinct classes of dynamical behavior. 


direction 

rate 

Model A 

Model B 

right 

k+ 

e^F/x (l-e2)/2 

e^Fjz ~ 1 - 2e 

left 

k- 

e-^Fjx ~ £^2 

e-^Fjz ^ £2 

up/down 

k° 

1/4 

1/Z~ £ 


TABLE I. Jump rates of the particle on the square lattice with 
periodic boundary conditions. Here X = 2(e^^^ and 

Z = -I- -I- 2, so that all rates = 1/4 for F = 0. 



FIG. 1. (Color online) Sketch of the two-lanes toy model, and 
its mean velocity [see 01 as a function of the force, for model 
A and model B. 


In order to have a quick taste of the qualitatively dis¬ 
similar mobilities that model A and B may give, let us 
introduce a minimal two-lanes model of particle diffusing 
in a meandering channel. This model includes two kind 
of states: C\ states where the force pushes the particle 
against a barrier (colored boxes in Fig[^ and C 2 states 
where the force pulls the particle away from a barrier 
(white boxes). One easily finds the density of states p 
and the mean velocity of particles in the steady state. 


p(Ci) = fc+ + fcO, 

p(C2) = k~ F k°, V =fc+p(C2) - k~p{Ci). (2) 


For models A and B the mean velocities are 

em _ e-^/2 


1 1 F 

VA = -tanh—, 


VB = 


(e^/2 -h e-P/2 -h 2)2 ^ 


( 3 ) 


respectively. The first is clearly increasing with F while 
vb cannot be monotonous in F as ub ~ e for F —> 00 
(Fig.[T] shows the plots of both functions). Hence, Model 
A does not yield NDM in this case, but Model B does. 

A similar fact pertains to scaling theories focused on 
the density of obstacles [3 HU]. For example, the theory 
presented in Ref. m does not yield NDM if constant 
transversal rates are plugged in their Eq. (2) for the av¬ 
erage velocity of the tracer, rather than those ~ e of 
model B. Indeed, with model A rates, the same formula 
becomes a monotonously increasing function of F and 
NDM disappears. 

In the following section we show a similar effect for 
a variant of the two-dimensional diffusion within slowly 
moving barriers, which was recently used PITU] as a main 
benchmark for testing ideas on NDM. 


III. SINGLE TAGGED PARTICLE WITHIN 
FLOATING BARRIERS. 

Model A was used [5] (with a slightly different nor¬ 
malization of rates) to show how NDM arises for par- 
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tides moving in crowded environments when the (non¬ 
overlapping) obstacles have low enough mobility 7 : given 
a configuration C with N obstacles and 1 tagged particle 
[see an example in Fig.j^a)], the bare escape rate from C 
is A(C) = Ai-|-A^ 7 , where Ai = + k~ +2k^ is the escape 

rate of the particle from any site within an empty grid 
(Ai = 1 for the rates summarized in Table |T]). We extract 
a waiting time r according to an exponential distribu¬ 
tion p{t) ~ exp[A(C)r]. After the extracted r, a move is 
attempted, picking up the tagged particle with probabil¬ 
ity 1/A or a random barrier with probability 7 /A. The 
tagged particle is then moved with a random selection of 
the direction that follows the rates previously exposed, 
while a barrier tries to jump with probability 1/4 to a 
random nearest neighbor. An update realizes the new 
configuration only if it is allowed, otherwise the config¬ 
uration is unchanged and another waiting time is added 
to the already spent r until a valid new configuration is 
achieved. Such modified Gillespie algorithm is very 
simple to implement and allows us to sample efficiently 
the dynamics of these models. Note that a time scale 
T* = 1/7 is associated to barriers’ motion, as opposed to 
a time scale equal to 1 for the particle motion. 

We have performed simulations within square lattices 
of side L filled with density pb = N/L'^ of floating barri¬ 
ers. An example of a typical configuration for pb = 0.2 
is shown in Fig. [^a). In all simulation we find that both 
model A and B display a similar NDM at sufficiently 
large values of F. 

However, it is not possible to switch with ease from 
one choice of microscopic jump rates to another. In¬ 
deed, even with the same crowding of obstacles, the sit¬ 
uation changes drastically if these cannot share corners 
[Fig.j^b)] and thus cannot join to form concave traps. In 
this case model A looses NDM while model B continues 
to display NDM (Fig. because the tracer cannot eas¬ 
ily move away from a situation as in Fig. j^b), where the 
rates of model B lead to an effective trapping of the parti¬ 
cle. By now this should be hardly surprising, as the con¬ 
straint of non-overlapping barrier’s corners leads to mor¬ 
phologies analogous to those exemplified by the meander¬ 
ing channel model of the previous section, where no con¬ 
cavity is present and model A finds no traps to generate 
NDM. In the limit of low pb the corner constraint should 
be irrelevant. Yet, theoretical approaches [3 uni were 
used to compute a nontrivial mobility also in the limit 
of low obstacle densities. Our results suggest that this 
procedure leads to an incomplete picture if only model B 
is considered. 


IV. MULTIPLE TAGGED PARTICLES WITHIN 
FLOATING BARRIERS. 


The crowding of an environment may result from the 
presence of different molecular species, and in general it 
is not straightforward to anticipate the mobility of par¬ 
ticles in complex conditions. Here we may think of the 
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FIG. 2. (Color online) (a) Configuration with density of ob¬ 
stacles pb = 0.2. The circle (red online) is the tagged particle 
and the squares are the barriers. The four obstacles marked 
with a cross are forming a trap for the tracer, (b) Same 
pb = 0 . 2 , when obstacles cannot have vertices overlapped. 



FIG. 3. (Color online) Mean velocity vs. force for models A 
and B with non-overlapping corners [case (b) of Fig. i, for 
L = 20, pb = 0.2, and two values of the obstacles’ mobility 
7 . Similar results are found for other small 7 ’s, in particular 
model A is quite insensitive to the value of 7 . 


supercrowded environment of the cell [30] which affects 
the diffusive properties of passive particles |3T]. Much 
less is known about the behaviour of active particles like 
molecular motors in such an environment |32j . 

As a first step toward more complex systems, in the 
simple model described in the previous section we may 
think the crowding not only determined by its barriers 
but also by a finite, possibly large number of tracers. 
The total density of occupied sites in the system is thus 
the sum of the density of barriers pb and the density of 
tracer particles Pp. The resulting model is, more than 
that of uni, a combination of the symmetric and asym¬ 
metric exclusion processes. The condition of a finite yet 
low density of tracers seems similar to the condition of 
low filling of bands of electrons, which leads to NDM in 
low temperature semiconductors unmi]. 

For model A it is not difficult to find conditions of low 
p = Pp + Pb where NDM is present as well as conditions 
with high p where NDM is not observed, see the exam¬ 
ples in Fig. [^ The latter state of normal mobility arises 
because the behavior of tracers emerges as an average of 
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FIG. 4. (Color online) (a) Configuration for model A with 
L — 20, Pp = 0.04, pb = 0.25, 7 = 10~^, and F = 8 pushing 
particles to the right. Tracers that may follow F by eventually 
never stepping back against it are drawn as empty circles 
(only one in this case), (b) The same for pp = 0.08. (c) 

Mean velocity vs. force (same parameters) for several particle 
densities. 

two typical conditions, one of tracers stuck in traps, and 
one of tracers free to move because all available traps are 
occupied by other tracers. For several values of pb we find 
that the contribution from running tracers determines a 
normal mobility for sufficiently high values of pp. This 
is expected to occur when there are more particles than 
the average number of traps [e.g. as in Fig|^b)]. If in¬ 
stead Pp is low enough, trapping is still the dominating 
behavior [see Fig^a)] 

This model therefore is an example where increased 
crowding leads to the disappearance of NDM. Again, this 
strengthens our point that the link between crowding and 
the mobility of particles is not straightforward. On the 
other hand, the numerical results may be correctly in¬ 
terpreted at least qualitatively in terms of saturation of 
trapping effects. Also in this case one could devise sim¬ 
pler systems where the mechanism of filling of all the 
traps is more easily detected. For example, in the two- 
lanes models where one lane is occupied by “hook” states 
that act as traps mil!, we have numerically verified that 


NDM is present only if the number of tagged particles is 
at most equal to the number of traps. 

V. DISCUSSION AND CONCLUSIONS. 

We have presented results from some simple models 
suggesting that the phenomenon of NDM in many cases 
can be ascribed to the presence of long living traps that 
catch the strongly pushed traced particles. There is in¬ 
stead no one to one correspondence between crowding 
and NDM. In particular, according to our results, it 
seems far fetched to draw conclusions from either sim¬ 
ulations or from mesoscopic theories if these rely on the 
specific choice of the jump rates. 

The choice of transversal jump rates normalized so that 
they decrease with the force is quite popular but it needs 
to be understood and justified. For example, starting 
from the concept of Brownian motion of a particle in a 
fluid, we do not see how the increasing force could lead to 
a weaker transversal motion, namely to a diffusion coef¬ 
ficient in the fluid that depends on a transversal F. Also 
for particles diffusing by jumping in a pattern of energy 
minima it is not general to find a decrease of transversal 
motion generated by an applied force. Think for exam¬ 
ple of a particle density within a minimum, at the micro¬ 
scopic level of description (while the mesoscopic level is 
the one of jump systems describing the energy minimum 
as a single state). If the applied force shifts the maximum 
of this density closer to the energetic saddle points that 
connect to nearby states, also in transverse directions, we 
would witness a case where the force would actually en¬ 
hance the transversal motion. It seems that such scenario 
is in principle as likely as that where the transversal mo¬ 
tion is depressed by the force. Hence, transversal jump 
rates that decrease exponentially with the applied force 
should be tested in parallel with other choices, keeping 
in mind that any specific system might be a realization 
of one of these choices. 

Our examples also show that focusing on the density 
of objects responsible for the crowding is not a stand¬ 
alone strategy. One needs first to determine whether the 
tagged mobile objects (e.g. particles or polymers) might 
become more stuck by the crowded environment when 
the pushing force is increased. Furthermore, it is also 
relevant to check if on average the traps are at least as 
many as the tagged particles. The study of NDM in more 
complex simulations than those considered so far (such 
as the condition outlined in the previous section where 
multiple tracers coexists) should help further clarifying 
these issues. Of course this is just part of a broader sce¬ 
nario, where kinetic constraints, jamming, glassy dynam¬ 
ics, alternating external forces, self-propulsion, etc. fur¬ 
nish many more mechanisms leading to negative mobility. 
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